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Введение. Рассмотрен частный случай синтеза одной и той 
же электромеханической системы управления методом 
максимума Понтрягина и методом синергетического 
синтеза. В качестве задачи была определена задача 
синтеза оптимальной по быстродействию 
электромеханической системы позиционирования, при 
этом момент сопротивления движению линейно зависел 
от выходной координаты системы. Этот подход к выбору 
задачи синтеза был обусловлен тем, что синтез 
оптимальных по быстродействию систем является одной 
из самых широко распространенных задач, которая 
решается при повышении эффективности действующих 
систем управления. 

Материалы и методы. Синтез оптимальной по 
быстродействию линейной системы управления на основе 
принципа максимума — широко распространенная задача 
в современной теории управления. Однако процедура 
синергетического синтеза такой формализации не имеет. 
Исходя из этого, в статье предложен подход, сближающий 
эти два метода, который, по мнению авторов, позволит 
повысить эффективность метода  синергетического 
синтеза, добавив в него некоторые особенности 
методологии синтеза оптимальных систем. 

Результаты исследования. В работе сформулированы два 
основных научных положения. Первое - применение 
принципа максимума для объекта класса двигатель 
постоянного тока при синтезе алгоритма 
позиционирования в условиях линейной нагрузки, 
функционально зависящей от угла поворота двигателя, 
позволяет оптимизировать систему по быстродействию. 
Второе - синтез системы управления на основе 
синергетического подхода позволяет получить систему, 
близкую к оптимальной (квазиоптимальную), но уже 
после модификации самого метода синергетического 
синтеза. Сформулирована гипотеза о возможной связи 
реализации процедуры 
состояния в методе 
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временем переключения управления. 
Обсуждение и заключения. Выполнен синтез управления 
методом максимума и методом Аналитического 
Конструирования Агрегированных Регуляторов АКАР. На 
основании сравнения эффективности применения методов 
выдвигается гипотеза о возможной совместимости 
исследуемых методов. 
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Введение. В 20-м столетии потребности техники, в частности космической, выдвинули серию задач, 
для решения которых была разработана новая теория, получившая название теории оптимального управления 
[1]. Одним из основных методов при синтезе оптимальных систем управления является разработанный в 50— 
60 гг. двадцатого века метод максимума, над которым работали советские математики Л. С. Понтрягин и его 
ученики [2]. Применение этого метода основано на формализации задачи синтеза с переходом к форме задачи 
Майера и последующем решении систем линейных или, в некоторых особых случаях, нелинейных 
дифференциальных уравнений [1]. Стоит отметить, что необходимость решения систем дифференциальных 
уравнений, а в нелинейном случае это не всегда возможно, во многом является ограничением, как самого 
метода максимума, так и всей концепции синтеза оптимальных систем управления. 

В последние десятилетия широкое распространение получил новый подход к синтезу систем 
управления, который в своей основе опирается на синергетическую концепцию анализа и синтеза систем. 
Используемый метод, автором которого является А. А. Колесников, получил название метода аналитического 
конструирования агрегированных регуляторов (АКАР) [4-5]. В основе предложенного подхода лежит 
концепция синтеза нелинейных обратных связей. Они обеспечивают асимптотическую устойчивость системы 
управления по отношению к требуемому движению аттрактора в пространстве состояний [4]. Отличием этого 
метода от методов синтеза оптимальных систем управления является то, что здесь отсутствует как критерий 
оптимизации системы управления, так и констатация оптимальности синтезированного управления. С точки 
зрения реализации процедуры синтеза, метод АКАР обладает несомненным преимуществом перед методами 
синтеза оптимальных систем, которое выражается в отсутствии ограничений на нелинейность системы 
дифференциальных уравнений [4—5]. 

С практической точки зрения синтез управления процессами или объектами должен обладать 
возможностью ответа на вопрос, а есть ли другое управление, которое обладает свойством превосходства над 
всеми другими. Так современное массовое производство, в целях обеспечения конкуренции на рынке, 
постоянно требует повышения эффективности. Одним из наиболее востребованных способов такого 
повышения эффективности является оптимизация процессов управления. Этот подход применим как к 
системам автоматизированной сборки оборудования [6-7], так и к металлорежущим системам, в частности к 
системам управления сверлением в металлорежущих станках [8-15]. Идея совмещения этих подходов к синтезу 
систем управления обладает, с точки зрения авторов, несомненной практической ценностью. С научной точки 
зрения актуальным является совмещение преимуществ метода АКАР, которые выражаются в возможности 
учета при синтезе сложной нелинейной динамики процессов в машиностроении [16—19], с четкой и ясной 
формализацией постановки задачи синтеза и оценки ее достижения в методе максимума [1—2]. 

1. Синтез базовой математической модели и постановка задачи исследования 

В современной экономике широкое распространение получило направление синтеза оптимальных по 
быстродействию систем, которое позволяет добиться существенного роста эффективности систем в целом [7— 
8]. Так, при решении задач автоматизации процессов сборки различного оборудования, часто возникает задача 
крепления различного вида частей друг к другу [6-7]. Здесь от скорости проведения этой операции зависит 
экономическая эффективность всего процесса производства. Эта же ситуация наблюдается в системах 
обработки металлов резанием на металлорежущих станках [9-19], в которых чем быстрее осуществляется 
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процесс обработки деталей, тем меньше издержки всего производственного процесса. Исходя из рассуждений, 
примем в качестве допущения, что требуется синтезировать оптимальную по быстродействию систему 
управления, которая с учетом схожести процессов автоматизированного крепления деталей и сверления 
глубоких отверстий может быть проиллюстрирована следующей схемой (рис. 1). 
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Рис. 1. Иллюстрационная схема 


На рис. 1 представлен двигатель постоянного тока, осуществляющий либо операцию крепления 
оборудования, либо подачу в сверлильном станке. В случае операции крепления — С суммарная 
характеристика жесткости сжимающегося при скручивании материала. Для случая сверления — С 
характеризует линейное нарастание сопротивления подаче при накоплении стружки в стружкоотводящих 
каналах сверла. 

В основе исполнительной подсистемы системы управления в обоих рассматриваемых случаях лежит 
электропривод, обеспечивающий преобразование электрической в механическую энергию движения 
исполнительных механизмов системы. Примем, что в рассматриваемом случае такое преобразование 
осуществляется электродвигателем постоянного тока с коллекторным управлением, который описывается 
следующей системой уравнений [20]: 


И-сю= ри + № 
ГИ , 
с = о. 
ПИ 
где И — напряжение, подаваемое на коллектор двигателя, т — ток, потребляемый двигателем; К, Е — 
параметры электрической части двигателя; Л — параметр, характеризующий инерционные свойства ротора 
двигателя, приведенный инерционный момент всех вращающихся масс; «о —[ частота вращения ротора 
двигателя; М.— внешний приложенный момент сопротивления; Си, се —— механическая и электрическая 
постоянные двигателя. В рассматриваемом случае момент представлен как функция угла поворота ротора 


двигателя, то есть М, = Со. С учетом этого (1) перепишем как: 


И. 
И 














: ао 
сСй=Л—+Са . (2) 
ИИ 
ао _ 
ИИ 
Выразив во втором уравнении значение тока и подставив полученное в первое уравнение, получим: 
1 [7 45 Е а® 1© КС 
О—-®= + + ®+ [о 
С. С. Ст Ч ССт Ч С. Ст ССт # (3 ) 
аа, 
О 
й 
Разрешив первое уравнение относительно старшей производной и переписав его вторым, получим: 
49, 
ж “ 
! 4 
4`®  Вао сс С ВС с’ ый 
5 = +—)о о+—- 0 
Ге Г 4 Гл У ЕЛ ЕЛ 
а® 
Обозначим переменные Ч=х,®Ф=Х,, = х, ; В качестве постоянных введем 
1 
сс, С ВС с 
=а,,,([——+—) =а.,, =а,,—” =Ь. Тогда система (4) примет вид: 
Г М Л М Ш 
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Или в матрично-векторном виде: 


х=Ах+ Ви, (6) 
от о 0 х 0 
где А=| 0 0 |, В= 0], х=|х, | и= 0 
—@: —@) @з Ь Хз и 


Примем в качестве привода, обеспечивающего движение рассматриваемой системы, электродвигатель 
постоянного тока, параметры которого обеспечивают следующие значения постоянных: аз!=4,65, аз›=4,6, азз=2. 


0 1 0 
Матрица А, с учетом этого, примет вид: А=| 0 0 1 |. Зная матрицу А, найдем собственные числа А 
—,65 —,6 —2 


матрицы А: А. =-1,2666 + 0,00001, ^› = 0,3667 + 1,88061, Аз = 0,3667 — 1,88061. 

Как видно из полученных собственных значений матрицы, система управления асимптотически 
устойчива по Ляпунову [21]. Более того, можно утверждать о колебательном характере процессов протекающих 
в системе, так как собственные значения не только содержат отрицательные вещественные части, но имеют 
ненулевые мнимые части. 

2. Синтез управления методом максимума Понтрягина 

Для начала сформулируем задачу синтеза оптимального управления в следующем виде: 

— используя принцип максимума для объекта, описываемого Системой (5), определить алгоритм 
оптимального уравнения, который обеспечивает перевод объекта из начального состояния х/(0)=50, х2(0)=0, 
х3(0)=0 в конечное состояние х/(Т)=0, х>(Т)=0, х(Т)=0 за минимальное время Т. При этом необходимо 
определить число переключений, моменты переключений и построить кривые управления и(1) и координат х/(1), 
х2(1), хз(и). 

На управляющее воздействие наложено ограничение |и<О„„. Параметры  системы:. 
а, =Ьа, =12, а. = 20,6 =1; Оак=440 В. 


Решение: Запишем гамильтониан: 


н-У чл. ® 





где лу; и } — координаты векторов \уи }. В дополнение к системе (5), составим систему уравнений для 


вспомогательных переменных \у,, где = 1, 2 ‚› используя следующую зависимость: 


2 
а = Уи, 9, о $=1,2. 





Или в раскрытом виде: 














И 








—- и т 
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С учетом (5), система уравнения (8) примет вид: 











а+ 
* =4,65$, 
а+ 
г =—Ф, +4,6$, (9) 
$ 
- =, +24, 
ГА 
Матрица собственных значений для случая (9): 
1,2666 0 0 
а 0 0,3667 + 1,8806 0 
0 0 0,3667 +1,8806 
Решение системы (9) для диагонализированного случая матрицы А%-Р примет вид: 
ф = Се'266ы 
т 
фр = С е(0,3667+1,88060! (10) 
2 2 ы 


р (0,3667-—1,88061)1 
$. =С.е 


С учетом матрицы У* ‚в исходном базисе решение примет вид: 
4, =-0,9474Се'?°%" + (0,5812 -0,33470С,е 69718806! + (0,5812 + 0,3347) Се 769718805 
м7 =-1 1893С е''2666! +0 6883С е(©;3667+1,88060! _0 6883С е(©.3667—1,88061: 
2 ы 1 В 2 ? Е ? 

ф, =-0,2581Се"?°°”" + (0,1812+0,2086С, ет + (0,1812 -—0,20861) Се 2967 18806 
где Су, С», Сз — постоянные интегрирования. 
Общее выражение, описывающее функцию Гамильтона: 

Н = +В (11) 

В выражении (11) важную, с точки зрения метода синтеза, роль играет член, включающий в себя 

управление: 
Н' =. 0 = (-С, 0,258 1е'°®°' + С, (0,1812+0,20860 + С, (0,1812—0, 2086 е®°°°71*8%° 9) . 

Для того, чтобы гамильтониан Н, вычисляемый по формуле (11), принимал максимальное 
положительное значение, слагаемое Н’ должно быть всегда положительным и наибольшим. Для этого алгоритм 
оптимального управления должен имеет вид: и(1)=с Иа», где 

с = 5421 (—0,2581Се'"°%°" + (0,1812-+0,2086е°7°°7 856 + ((0,1812—0,2086е®.°5°7 18805’) 


Оптимальное управление и(1) является кусочно-постоянной функцией, принимающей значения = Оша», И 





имеет не более двух интервалов постоянства, так как нелинейная функция 

—0,2581С е”?°°°' + (0,1812-+0,20867е°567 8060: +. ((),1812—0,20861е®'2567 88060: , 

меняет знак не более одного раза. При этом возможная смена знака происходит с плюса на минус, то есть для 
выполнения принципа максимума требуется сначала подавать на двигатель Ч= +0», а затем Ч= —Ошах. 
Проверим эти рассуждения, построив полученную функциональную зависимость 4’, в пакете МаЧаьБ (рис. 2). 








Рис. 2. График изменения Ч’, 
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Как видно из рис. 2, 4, действительно меняет знак с (+) на (-) только один раз. С учетом этих 


рассуждений определим алгоритм оптимального уравнения, который обеспечивает перевод объекта из 
начального состояния х!(0)=0, х›(0)=0, хз(0)=0 в конечное состояние х!(Т)=50, х›(Т)=0, хз(Т)=0 за минимальное 
время Т. Здесь отметим, что в современных условиях нет необходимости получать аналитическое решение 
исходной системы уравнений. Можно, используя имеющиеся возможности современных пакетов прикладных 
программ типа МаЧаЬ, достаточно просто и наглядно получить численное решение для рассматриваемого 
случая. 

Результаты численного моделирования для расчетного оптимального по быстродействию алгоритма 
управления представлены на рис. 3. Параметры модели подобраны таким образом, чтобы она выполняла 
требуемое по правой границе граничное условие. 


Хь, рад 





0 0.5 1 15) ъ сек 


о, В а) 
500 
0 
т, сек 
-500 
0 0.5 1 1.5 
ь) 


Рис. 3. Графики выходной координаты (а) и переключения управления при 1 = 0,5 (Ъ) 


Как видно из рис. 3, время переключения управления с (+0 „„х) на (-Чвх) выбрано как 1 =0,5 секунды, 
а общее время регулирования составило и Т=1,3 секунды. При этом система управления выходит на 
требуемый нам уровень выходного сигнала, то есть достигается правая граница краевых условий. 

3. Синтез системы управления методом АКАР 

Для синтеза системы управления методом АКАР можно не осуществлять переход к абстрактному 
случаю пространства состояния (см. (5)), а удобно использовать исходную систему (2). Однако требуется 
преобразовать ее к следующему виду: 


ое (12) 


ь С С с, К 1 
где переменные а&=х,®Ф=х,,Ё=х,; в качестве постоянных введены > аз, а а Чт Ь.С 


7 21? Г, 32, Г, 
учетом определенных ранее значений, получим: а. =1,515, а, = 2,325, а =1, 5,аз =2,6 =1 и система (12) 











примет вид: 
ле 
2 
" = 2,325, +1,515х, , (13) 
ГА 
се, 
ГИ . 
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Матрица А, с учетом этого, примет вид: 


0 1 0 
А=| -2,325 0 1,515 
0 —15 —2 


Зная матрицу А, найдем собственные числа /). матрицы А: 1 = —1,2674 + 0,00001; ^› = 0,3663 + 1,88011; 
Лз = -0,3663 -— 1,88011. Как видно из этих результатов, речь идет именно о рассматриваемом при синтезе 
методом максимума случае системы управления. 

Координата, характеризующая угол поворота ротора (х1), является выходной координатой системы. 
Поэтому для формирования требований к желаемому поведению системы в пространстве состояний введем 
макропеременную порядка 4, =х,—х„ => 0, где х„ — заданное требуемое значение координаты х,. Значение 
скорости изменения угла должно обеспечивать стремление угла поворота ротора к заданному выше значению. 
Для этого введем еще одну макропеременную Ч =х, -ф.(х)=0, где ф,(х) — некоторая функция, 
описывающая стремление х, >ф,(х) в стационарном состоянии. Дальше вводим макропеременную порядка 


по координате -е =х, -ф.(х,х,) >0, где ф.(х,х,) — некоторая функция, описывающая связь между 


координатами в стационарном состоянии системы и, следовательно, х, => ф.(х,^х,). 
Для вновь введенных макропеременных Ч,, 4,, 4, потребуем асимптотически устойчивый закон 


изменения, то есть выполнение следующей системы дифференциальных уравнений: 








а+ 
Ти +ч =0 
а+ 
Т— 2+4, =0, (14) 
$ 
Т— +, =0 
где Т, Т,, Т. —Щ Постоянные времени, определяющие скорость стремления введенных макропеременных к 


нулю, либо, иными словами, скорость стягивания произвольной траектории движения системы в пространстве 
состояния к требуемому аттрактору. 

Введенная в выражении (14) система уравнений расширила пространство состояний системы с 3-го до 
6-го, что является одним из основных положений метода АКАР. Дальнейший синтез системы управления 
сводится к поэтапному процессу декомпозиции (сжатия) этого пространства до начального уровня [1-2]. 
Начнем процесс декомпозиции с подстановки в первое уравнение системы (14) вместо макропеременной ЧФ ее 
значения в координатах управляемого процесса. В результате получим: 


пу жи) 0. (15) 
С учетом (13), уравнение (15) примет вид: 
Тх-+@-х)=0. (16) 
Учитывая следующий шаг алгоритма синтеза, в результате которого приняли, что х = ф,(х,) ‚ найдем 


значение ф,(х), как: 


(хи) 
ф-(х)=- . т о (17) 
1 
Учитывая это, макропеременная $, примет ВИД: 
И с Ре. = (18) 


1 
Справедливость выражения (17), с точки зрения целей управления, подтверждает тот факт, что в стационарном 
состоянии значение х, > ф,(х,) > 0. Тогда, с учетом (18), второе уравнение системы (14), примет следующий 


вид: 


Ах аф.ь(х 
О, -9.(н)-0. (19) 


Или, с учетом значения ф,(х): 
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На следующем шаге алгоритма синтеза, пользуясь теми же рассуждениями, что и раньше, и тем 
фактом, что х, > ф.(х,х,) , определим из (20), с учетом (13), значение ф,(х,,х,) в координатах управляемого 
процесса: 


м 
3 1 д 1 2 


1 2:1 
Получив значение ф-(х,х,) в координатах управляемого процесса, можно определить значение 





о. 21) 


макропеременной ЧФ. в координатах состояния системы и раскрыть третье уравнение системы (14). 


Ш. а$ф.(х,х 
+, 4,50. (22) 
После подстановки в (22) полученных ранее значений ф.(х,,х,) из (21) и использования системы (13) 


т 


на этом шаге декомпозиции пространства состояния системы управления, определим значение управления в 
координатах управляемого процесса: 
(Т.+Т, +ТИ,)_ Ч, +ТЛ, +1) 
3 
ел =: Хо) 
"РТ ТТЛ 

Выражение (23) и определяет асимптотически устойчивое управляющее воздействие на систему 
управления, описываемую уравнениями (13). 

Результаты моделирования системы дифференциальных уравнений (13) с требуемым значением угла 





И =1,5х, +2х, +1,535х, +1,5х, 


(+Т,-Т,) р 


<> 0,66х, 





0,66 


поворота ротора двигателя х=50 рад и значениями введенных постоянных времени Т =0,1, Т, =0,2, 


Т, =0,3 представлены на рис. 4. 
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Рис. 4. Результаты моделирования с управлением при Т; = 0,1, Т, =0,2, Т, =0,3 : 





переходный процесс по выходной координате (а); график изменения управления (Ъ) 


Как видно из рис. 4, в сравнении с рис. 3, структура управления в обоих случаях одинакова, то есть 
речь идет о том, что в начале подается положительное по знаку управление, а затем его знак меняется на 
противоположный. Время регулирования резко увеличивается до значения близкого к 2,5 секунде, но при этом 
максимальное управление достигает значений больших 5 000 В. Это недопустимо с точки зрения метода 
максимума. Для ограничения максимального значения управления введем в контур управления нелинейное 
звено типа «насыщения» с таким же порогом, как и в случае максимума. Результаты моделирования системы с 
учетом этого порога и с подобранными значениями постоянных 1, Т,, Т, представлены на рис. 5. 
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Рис. 5. Результаты моделирования с управлением при Т, = 0,15, Т, =0,16, Т, =0,17: 


переходный процесс по выходной координате (а); график изменения управления (Ъ) 


Как видно из рис. 5, введение нелинейного ограничения на управление типа «насыщения» позволяет 
получить требуемый с точки зрения допустимости управления результат. При этом фактически соблюдается 
определенная методом максимума структура управления, но управление не является оптимальным по 
быстродействию. Здесь необходимо отметить, что такая модификация метода АКАР может приводить и к 
потере устойчивости системы (рис. 6). 
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Рис. 6. Результаты моделирования с управлением при 1 =0,015, Т, = 0,016, Т., =0,017: 


переходный процесс по выходной координате (а); график изменения управления (Ъ) 
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Однако, несмотря на результат, представленный на рис. 6, дополнение системы управления полученной 
методом АКАР, ограничениями на амплитуду управляющих сигналов существенно сближает ее с решениями, 
получаемыми при реализации процедуры синтеза методом максимума. 

Обсуждение и заключения. Таким образом, для базовой модели системы управления (2) было 
синтезировано оптимальное по быстродействию управление методом максимума. Для этого же случая, в 
третьей части работы с этими же параметрами системы и граничными условиями было синтезировано 
управление методом АКАР. Оба эти управления внешне несхожи. В случае метода максимума — это 
программное кусочно-постоянное управление, где важен оптимальный момент времени переключения знака 
управляющего сигнала. В случае АКАР — это непрерывная гладкая функциональная зависимость управления 
от значений координат состояния системы. Однако, при их реализации видна структурная схожесть, 
выражаемая в смене знака управления. После модификации метода АКАР структурная схожесть стала еще 
более видна. Отметим, что добиться параметрического подобия реализуемых управлений можно за счет 
подбора вводимых при синтезе методом АКАР постоянных времени. Это позволяет авторам сформулировать 
следующую гипотезу. 

— Подбор вводимых в процедуре синтеза методом АКАР значений постоянных времени 71,,Т,,Т, ‚в 
случае модификации полученного управления способом предложенным в статье, позволяет добиться от 
синергетического управления свойства оптимальности по быстродействию. Авторам не удалось получить 
общее доказательство этой интуитивно понятной гипотезы. Однако такая задача не стояла в рамках этой 
работы. В результате численных экспериментов с вариацией Т,Т,,Т, удалось добиться квазиоптимальности 


управления, полученного методом АКАР. 
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